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Abstract 


The interplanetary shock speed and surface normal for the July 6, 

1966 shock in the vicinity of the earth have been accurately estimated 
and shown to be in serious disagreement with previously estimated values. The 
shock speed was 57CH20 km/sec and the normal had a direction given by 
©se = -3 8° +7° and = l66°+7°. These values yielded the most 
consistent picture of the event where the following constraints were 
put on the estimation: (l) the MHD shock conservation equations (excluding 
the normal momentum and energy equations) had to be satisfied and (2) the 
time -of- sighting restriction for the shock passage between Explorers 
28 and 33 * both in interplanetary space, had to be fulfilled. Vela 3A, 
which was in the magnetosheath at the time of its shock sighting, observed 
the shock 3.0 +2.1 min. later than might be expected from the simple 
model of a shock front remaining plane and constant in both speed and 
(normal) direction after interaction with the bow shock. The lower 
end of this range (1 min.) is not inconsistent with the delay time 
obtained by Leaping and Chao (1971) in studying the January 11, 1968 
traveling shock-bow shock interaction. Hence, the traveling shock probably 
experienced a rapid slowing-down and significant distortion from a plane 
shape in the sunward region 0 f the earth* s magnet os heath. 



Introduction 
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The actual shapes of shock surfaces in interplanetary space have an 
obvious bearing upon models of the propagation of shocks caused by 
solar events. These models in turn influence explanations of the solar 
sourcs cf such shocks and yield knowledge of the detailed characteristics 
of the interplanetary medium. The July 8, 1966 shock /'Van Allen. 1966; 
Van Allen and Ness , I967, Ness and Taylor , 1968; and Lazarus and B insack , 
1968) in the vicinity of the earth has been estimated by Greens tadt 
et al. (1970) to have an interplanetary normal which is inclined 
65-70° below the ecliptic plane at solar ecliptic longitude 165°. They 
made this estimate using a three -spacecraft kinematic method. This 
direction is in the correct octant for a shock induced by a solar flare 
located at heliographic coordinates 34°N 5 48°W and occurring at 0022 UT 
on July 7. If the severe inclination of this estimated normal is 
correct, it would provide important quantitative information regarding 
the possible shapes of such shock surfaces in interplanetary space 
(e.g. Hirshberg. 1968; Hundhausen . 1Q71 : and Schatten and Schatten . 1971). 

It is the principal purpose of this study to obtain two estimates 
of the July 8 interplanetary shock normal in the vicinity of the earth 
by using two semi -independent methods, each considered to be of 
approximately equal. accuracy. 



Methods 


The estimated normals from these two methods must be, on the 
one hand, consistent, within the error due to fluctuations in the data, 
with a least-squares best fit of interplanetary magnetic field and plasma 
data (from Explorer 33) to a six-equation subset of the Rankine-Hugoniot 
equations. This method described by Lenning and Argent! ero (1971) 
will be referred to as Method 1. On the other hard, these estimates 
must be consistent with the times of observation of the shock by two 
spacecraft in interplanetary space; this requirement is the essence of 
Method 2. 

We now describe Method 2 in detail. As a by-product of Method 1 

-4 -4 -4 

the shock speed V s and the magnetic field difference AB(s Bg-B 1 ) across 
the shock are obtained. 

The components of the best-fit plasma and field quantities (including 
densities as one-component vectors) can be examined from the point of 
view of their estimated errors normalized by the magnitudes of their 
associated vectors. The source of the errors is primarily the natural 
fluctuations in the original data over the analysis interval as measured 
by RMS deviations. When these normalized errors are compared against 
each other, those for the field quantities are usually distinctly smaller 
than those for the plasma quantities; this is in fact the case for the 
July 8 shock (Table l). Also field differences across the shock are 
expected to have, in general, smaller errors than the fields themselves 

«“4 

because any original measurement -bias would tend to cancel. Hence, AB is 

probably the most accurately determined quantity from the best-fit analysi 

— » — * — * 

Then for known AB, V s , R and T (where H is the displacement vector 
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between any two interplanetary spacecraft and r the time delay between 

shock sightings for the two spacecraft) the shock normal is estimated 

by the following scheme : 

-» — » 

AB . n - 0 ( 1 ) 

-• -* 

R . n =TV S ( 2 ) 

M = 1 > (3) 


where the single Rankine-Hugoniot equation (1) is used and where constant 

-4 

n and V over the region of interest is assumed. For the second 
s 

spacecraft Explorer 28 (IMP 3) was used. Plasma data was not available 
from this spacecraft so Method 1 could not be applied in its case, 

— ► -4 

Since |r|/|t| ( £l6,000 km/sec) for Explorers 28 and 33 is so much larger 

than characteristic interplanetary shock speeds, n depends only very 

weakly on the value of V_ through the constraint imposed by equation 

(2). Hence, Methods 1 and 2 are coupled essentially only byAB, believed 

to be the most accurately determined quantity in the analysis. If the 

two methods give normals which agree within a reasonably small error 

cone angle, then the necessary self-consistency is accomplished. 

Vela 3A, which was used by Greenstadt et al. in their estimation 

of the interplanetary normal together with the above-mentioned two 

spacecraft, was located in the magnetosheath at the time of its shock 

sighting and therefore, we believe, disqualified to supply data for 

Method 2. It is not sufficient to say that Vela 3A was ''close to” 

the boundary of the bow shock and therefore almost in interplanetary 

space, because the traveling shock will slow down severely just after 

its encounter with the bow shock giving a large relative time-delay 

error ( Dryer et al. « 19&7 a^cl Lapping and Chao 9 1971) 5 this implies a 

large normal direction error in Method 2 • Also , as pointed out by 

Bur lag a (1971) , the three spacecraft were all close to the ecliptic 

plane, and therefore the component of the normal perpendicular to the 
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ecliptic plane is the one least accurately determined by the Greenstadt 
et al. three-spacecraft method. In this study we use the solar ecliptic 
coordinate system for representing both physical quantities and 

A 

spacecraft positions. This system is defined such that R is in the ecliptic 

A 

plane radially away from the sun, N is normal to the ecliptic and 
"Northward”, and T = N x R, which is tangential to the earth 1 s velocity 
about the sun. Using Method 2 we notice that even though R (=&1.9,96.6,3.4)R e ) 
is nearly in the ecliptic plane, A B (=( -1.9,7. 6, -4. 8)r ^calculated via the 
least squares best-fit Method 1, is about 31° below it. 

In applying Method 1 it turned out to be important that only a 
rather short interval of data, taken near the shock, be used. There is 
always a certain measure of unavoidable arbitrariness in choosing the 
proper interval to be used in a least-squares shock analysis, but this 
arbitrariness is significantly decreased by the evaluation of an 
associated quality index. This index is defined as the square root of the 
ratio of the total number of points of all shock parameters used in the 
analysis to the standard sigma-weighted least -squares loss function at 
convergence ( Lending and Argentiero , . 1971) . This index is usually greater 
than unity for characteristic interplanetary shocks provided reasonable 
sigma-weights are employed. For the July 8 shock the short analysis 
interval yielded an index of 1.2 whereas longer intervals gave distinctly 
lower values, usually below unity. After some necessary trial jand error, 
the analysis (short) interval chosen as most reasonable was; 
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5 points (6.5 min.) before shock 

6 points (7.1 min.) after shock. 


Magnetic field .data near the shock 


Plasma data near the shock Average of first 2 points (a« 

22.2 and 16.7 min., respectively) 

( before shock, excluding a spurious 
point occurring just before the 
shock (including the spurious 
point caused a drop in the quality 
index to the unacceptably low value 
of 0.86). 

First single point (10.5 min.) 
after shock. 

The magnetic field data indicated significant directional changes about 
seven minutes before and after the shock jump as listed above; this, in 
the final analysis, determined the field data interval. The plasma data 
that was used, guided by the quality index, turned out to be simply 
the closest reasonable values to the shock jump. 

Since there were no plasma data within the analysis intervals 
determined by the magnetic field data (except for the single spurious 
point just before the jump) , use of such data requires justification 
that the plasma points that were available, were still characteristic 
of the "true” values, before and after the shock, within estimated 
errors. As Table 1 shows, the plasma data was weighted in the analysis 
by rather large sigma-weights, which were based on RMS deviations in the 
data over +30 mins from the shock jump. The magnetic field data, in 
turn, was weighted with relate vely smaller sigma -weights baaed on the +7 
mins intervals from the jump. The (9 % certainty) error cone angle 
which is estimated for the resulting normal, and given below, will depend 
on all of these sigma- weights, and it will be an optimum estimate based 
in part on the relatively large estimates of the plasma sigma-weights. 

The important assumption made here then is that the sigma-weights put 
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on the plasma quantities are characteristic of the fluctuations close 
to the shock, i.e., over the interval for which there is no plasma 
data (at +11 mins), as they are for the +30 mins interval. 

The best-fit analyses was repeated ten times using the same (short 
interval) input data in each case except for a change in each one of the 
five plasma parameters, which was successively replaced by either its 
extreme maximum or minimum value based on the average +<7. In no case 
did the normal deviate by more than 3° from the best -fit normal 
based strictly on the average plasma values given in Table 1 as input. 

Also in no case did a field component change by even as much as 0,17. 

Even though this indicates, in this case, a weak dependence of the 
resulting normal on the plasma parameters, the influence of these 
parameters is still significant a,s the small error cone angle (11°) and 
high quality index indicate. 

Table- 1 gives the best-fit results of Method 1. The 
best-fit magnetic fields, before and after the shock, expressed in 

solar ecliptic 0 and 0 (B = (F , 0 gE , 0 gE )) were 
— ♦ 

= (22. 6r, -6.4°, 322.1°) and 

■■f a. 

B 2 = (20.77, -17.4 , 301.4°), respectively. These may be 
compared with those of Van Allen and Ness (1967) which were 

B^V.A.N) = (12.27, -10°, 310°) and 
— * 

B 2 (V.A.N) = (20.87, -l8 , 300°), where these estimates were 
derived from straightforward 13 min. averages. From the point of view of 
accurately estimating shock normals from magnetic field data a Lignificant 
difference occurs only for the pre-shock parameters 0 S e; otherwise the 
agreement is rather good. The difference is principally due to the 

difference in the length of the analysis intervals used and partly 
on the influence of the plasma data on the present estimates. 
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TABLE 1 

July 8, 1966 Explorer 33 Shock Parameters 


Parameter Average VaJLue 


b xr 

(y) 

-8.3 

b it 


9.2 

b in 


-1.2 

®2r 


-lo*4 

B 2 T 


17.0 

B 2 N 


-6.3 

W * 
W R 

(km/sec. ) 

69.6 

W T 


-2.9 

WN 


-38.4 

N 1 

(#/cm 3 ) 

5.3 

N 


8.0 

2 



°R 


0.83 

n T 


-o.i4 

“n 


-0.54 


Estimated Sigma Best Fit Value 


0.70 

-8.4 

0.70 

9.3 

0.75 

-1.4 

0.83 

-10.3 

0.57 

16.9 

0.75 

-6.2 

30.0 

55. C 

25.0 

20.7 

20.0 

-56.6 

0.49 

5.1 

0.78 

8.4 


0.85 


-0.13 


-0.52 


W is the difference velocity 



across the shock. 
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Tne two-point average of the preshock plasma velocity in an 

«■»> 

inertial coordinate system at r3st with respect to the sun was <V^> = 

(419, -3, -33) km/sec. Using this and the best fit plasma and normal 
values from Table 1 in the mass conservation equation yields an estimate 
for the shock speed V s = 570420 km/sec in the same inertial system. Then the 
component of the shock velocity in the ecliptic plane, based on Method 
l*s normal is, 490+17 km/sec, which is in close agreement with the 
estimate of this quantity bv Sugiura et al. (1968); using a different 
approach they arrived at a value of 500 km/sec. The error is due 
principally to the uncertainty in the magnitude of <V|> (+15 km/sec). 

Since Explorer 28 recorded the shock onset time and uncertainty interval 
be 2105:54+24 IJT ( Tavloy . 1969) and since the Explorer 33 onset time 
was 2105:40.5 (+3.5/-4.5) UT ( Greenstadt et aj . . 1970), the delay time 
^(—tgg-t^) between the two spacecraft ranged between -l4 and +42 
sec. Hence, tV s in (2) is somewhere between -1.3 Rp and +3*9 R-g, where 
the worst-case value of V g (590 km/sec) was used. We employ Method 2 
twice, once each for rV g equal to these two extreme values. Table 2 
lists the results of these calculations as well as the results of the 
estimates of other investigators for 
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TABLE 2 


Estimated Normals for the July 8 , 1966 Shock 


Investigator(s) Method n = ( n p> %>> n w) 9 rf . 9 $ re 

Van Allen and Ness Coplanarity Theorem 

( 1967 ) using Explorer 33 

average magnetic 
fields , GSFC experimei 

O. 89 , +0.03, -0.45 -2 7° +5° , 182° +5 ° 

Ness and Taylor Magnetic Field data 

(1968 from each spacecraft 

(Explorer 28 , Explorer 
33 respectively) using 
=0 and delay-time; 
coordination. 

I 0 . 90 , -0.33, + 0.28 16 °, 160 ° 

O. 65 , - 0 . 20 , -0.73 -47 , 163° 

i 

Taylor ( 1969 ) Coplanarity Theorem j 

using Explorer 28 J 

average magnetic j 

fields, then assum- j 

ing Q gE = 0 1 

1 O. 85 , -0.34, +0.41 24°, 158 ° 

j0.93, -0.37, 0.00 Q o 158 ° 

1 

■ 

j 

Greenstadt et al, Three spacecraft 

( 1970 ) method 

j 

,0.33, -0.09, - 0.94 - 70 °, 165° 

Greenstadt et si. Copianarity Theorem 

(iqjQ) using average magneti 

* ' fields. Ames experi- 

merit. Explorer 33 

Burlaga (1971) Greenstadt et al. AB 

field values, their 
shock surface inter- , 
section with the ; 

ecliptic, and using ■ 
AB*n = 0. 

Lepping Method 1 

10.78, -0.17, - 0.60 -37°, 168 ° 

3 

. 

0.74, - 0 . 21 , -0.64 -40°, 164° 

0 . 85 , - 0 . 13 , - 0.52 -31°, 171° 

Lepping Method 2 /tV =-1.3 R E 0.72, -0.23, -0.65 -4l°, 162° 

(tv|= 3.9 R e 0.77, - 0 . 19 , - 0.61 - 38 °, 166 ° 


Lepping 
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comparisons. The two normals obtained from Method 2 differ by less than 

4°, Obviously the small difference between the results of Methods 1 

and 2 is due to the inaccuracy in &B, The (95 $ certainty) error cone 

o 

angle associated with Method 1 was 11 . The angles between the normal 

o 

obtained from Method 1 and the two normals from Method 2 were 11 
(tV, = -1.3 R E ) and 5° (rV s = 3*9 R E ), respectively. Therefore, the 
two methods agree to within 11° with approximately a 95$ certainty. 

The averages of 0 gE and for these three normals are 

< e sE > = -37.7° and <0 SE > - 166.3°, 

which fall between the angles obtained by Greenstadt et al. , using only 

±he_ coplanar i t.v the orem and straightforward averag e magnetic field 

values . and the angles obtained by Burlaga. Tavlor 1 s (1969) assumption 

that © SE = 0° was not applicable in this case— also there was a very 

large uncertainty in the Explorer 28 Bg^ component in any case. 

Since the error cone angle associated with Method 2 is estimated 

-* o 

under worse case assumptions concerning the A B error to be also 11 , 

the average normal given by<n >=(<© SE > , <0 SE > ) has an estimated 

o ✓ 

error angle of about 7 . The Alfven Mach numbers, based on the best-fit 
magnetic field components along the shock normal <n >, were 3.0 + l.o 
(preshock) and 2.7 + 1.1 (postshock). These are somewhat lower than 
typical values for interplanetary shocks (**5) Hundhausen, 1970). 
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Further Results and Discussion 


Using (2) again with R (*AR y ) now representing the displacement 
vector between Explorer 33 and vela 3A, t v the related time delay, and 
V s equal to 570 km/sec, enables us to determine the time Vela would have 
seen the shock had the spacecraft been in interplanetary space 
during the shock transit. Where £R y ( = R 33 -R v ) was ( 39s 6l, 13) R £ 

we find that t v is 130 sec. Then the "expected" time of sighting at 
Vela, "t v " (= ^33"%) * s 2103:30. The actual Vela sighting took 
place at t v =2106:32 +1 according to Greenstadt et al. The difference 
between actual and "expected" times, 3*0 mins. + A #r v , where AT y is the 
uncertainty in the estimate of T v , represents the delay experienced by 
the traveling shock (T.S) due to the T.S.-bow shock interaction. The 
quantity AT V is 2.1 mins, for a 7° uncertainty in <n >, and the lower 
end of the range 3.0 +2.1 min is a reasonable value for such an interaction 
delay time (Leaping and Chao . 1971). The normal n y =(0.70, -0.25, -0.67) 
(0 SE = -42°, (zLm = l6o°) would, in fact, give such a result. The other 


end of the range, 5.1 min., is obtained from n y = (0.84, -0.11, -0.53) 
(0 qE = - 32°, - 173°). If the interaction delay time is indeed 


close to one minute, the: n v might represent a refinement of <n >(+7 ) 
but such a conclusion is too speculative without detailed knowledge of 
at least where Vela-3A was with respect to the bow shock. Notice 


however, that xiy lies very close to the normal derived by Method 2 for 

TV = -1.3 @*2° difference) . 
s 

The instantaneous shock speed in the vicinity of the earth was 
estimated by Van Allen and Ness (1967) to be 890+40 km/sec which is 
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probably 1-1/2 times too large. Our estimate of 570+20 tan/sec 

satisfies a consistent picture as described above, and furthermore 

is very mu *.h in line with all of the most accurately determined shock 

speeds listed by Hundhausen (1970). The Van Allen and Ness estimate of 

V was based on the assumptions that (a) the shock front was orthogonal 
s 

A 

to the radius vector to the sun (i.e., the R-direction at the earth’s 
location), (b) the onset time of the terrestrial sudden commencement 
(2102.2 UT) corresponded to the passage of the shock front past the 
center of the earth, and (c) the nearby moon had no significant effect. 
Assumption (c) is probably a good one and was implicitly used in this 
work. Assumption (a) is not acceptable since such calculations using 
time delays between ’'spacecraft" (earth and Explorer 33) are usually 
sensitive to the true direction of n, and R and < n>were 4o° from each 
other in this case. However, if a more accurately determined norm a l 
such as < n> = (-37*7 , ,166*3 ) had been used, too low a value of V 

s 

(;s? 380 km/sec) would have resulted. Using n v yields a still 

lower value (210 km/sec). This inconsistency is not surprising since 
assumption (b) is a weak one in light of the fact that the signals 

arriving earliest at the earth, which are initiated at the outer boundary 

of the magnetosphere R=Rm> propagate downward with characteristic 

speeds usually higher than the interplanetary shock speed and take about 

1-1/2 min. for R M = 10 R £ (Sugiura, 1965)5 these start as hydromagnetic 

waves and then degenerate into ordinary electromagnetic c missions below 

the ionosphere. Adding to this the fact that the traveling shock will 

slow down as it interacts with the bow shock makes assumption (b) 


difficult to accept. 
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The average shock speed from the sun to the earth, for the solar flare 
pccurring at 00 22 T on July 7, was determined by Van Allen and Ness to be 
950 km/sec with about a 1 % uncertainty. As is commonly the case this is 

considerably higher than the instantaneous speed at 1 AU ( Hundhausen . I970). 

Ness and Taylor (1968), using their estimate of the shock normal 

(0 SE = 16 , 0 SE = l60°) derived from the magnetic field data of Explorer 
28 and the near simultaneity of the shock onset at Explorers 28 and 33, 
calculated an upper limit on V s . Their estimate of this limit was based 
on the assumption that the onset time of the SSC at earth occurred 
(shortly) after the arrival of the shock at the bow shock. For the time 
delay (210 sec) and associated displacement vector in our equation (2) 
they used those quantities for Explorer 28 and the shock- shock interaction 
point. They point out that the propagation velocities through the 
magnetosphere and magnetosheath are significantly different from those 
in interplanetary space, and imply that a limit on V s is all that may 
be obtained without use of plasma data. The upper bound turned out to be 
710 +50 km/sec. If a more accurate normal such as<rf>is used in this 
model, essentially the same value (720 km/sec) is obtained. However, a 
normal mid-way between <A>and the Ness and Taylor normal yields a bound 
of 110 km/sec in this model, making their result somewhat fortuitous. 

For an accurate normal this sort of moclel should yield a reasonable but 
conservative upper bound on V_ for the type of magnitude of the time 
delay encountered here. 

Instead of estimating V s as these authors do, we "reverse" the 
calculation and proceed as we did for Vela where now the earth replaces 
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Vela in the calculation. The resulting time is an "expected” time 
that the earth would have seen the shock had its environment not 
"interfered" with the traveling shock. Using V g = 570 km/sec, <n>, and 
R e -33 = (29* 3 » 59*9? -2.3) R^ the effective delay time is seen to be 
T e _22 =2.3 min. This, of course, assumes that V g and n are constant 
between the earth and Explorer 33* Although crude this is a more realistic 
assumption than assuming these quantities constant in the magnetosheath 
in front of the earth where most of the slowing down of the travel ing 
shock takes place. Hence, the "expected" onset time at earth is 
2103.4 UT (because the onset time at Explorer 33 was 2105.7). Since the 
SSC at earth occurred at 2102.2, the time difference (t D ) of +1.2 min. 
between "expected" and actual onset times represents a net speed-up 
of the overall signal, regardless of its nature, as received at the 
earth. Sugiura et al . (1968) have shown, in fact, for the July 8 
event that the magnetospheric propagation of the magnetic field toward 
the tail is indeed faster than the propagation of the interplanetary 
shock just outside the bow shock, even if it is assumed that V a = 700 

o 

km/sec. They used the magnetic field data of OGO 3 which was located 
at (10.1, -10. 3 > 11.6) Rg. They also point out that a 60 sec. transit 
time from the magnetopause to the earth along the earliest-signal-path 
is not unreasonable for this event. 

Using this 60 sec. transit time let us estimate tp from theoretical 
considerations. It is clear that, according to our model, 


%) ~ ^SH " + 


(4) 
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where is the time the shock would travel from the bow shock to the 

earth (* 15 R^) at 570 km/sec unimpeded, t^g is the minimum hydromagnetic 

wave travel time of 60 sec. through the magnetosphere («12 R^), and t 

JKi M 

is the magnetos heath shock travel time where the shock velocity will be 
assumed equal to ^400 km/sec in this region (**3 R^) , Hence, t gH s2.8 
min., t^^0.8 min., and therefore 1.0 min., which compares 
reasonably well with t^ =1.2 min. above. Again there is considerable 
uncertainty associated with tp -observational and t^ -theoretical, but 
either value cautions against modeling the true t^ as zero. 
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Concluding Remarks 

The most accurate!, y estimated interplanetary normal for the July 8 
shock in the vicinity of the earth is indeed rather severely inclined 
with respect to the ecliptic plane. But an inclination of 0 gE s-38° 
is more consistent with a reasonable overall picture then the value of 

o o 

9 = -70 as obtained by Greenstadt et al. The -70 inclination would 

bb 

not allow a satisfactory fitting of the shock data to the conservation 
equations because this direction for the normal lies outside the 11° 

•4 

error cone whose axis is along 4n^ (or n from Method l) • The use of 
plasma data from Explorer 33 enabled us to estimate 

the shock speed to be 570+20 km/sec, significantly lower than the Ness 
and Taylor upper limit of 710 +50 km/sec. The onset time of signals 
at Vela and the earth aided in more fully understanding the shock 
interaction event: apparently the traveling shock experienced a rapid 
slowing-down and distortion from a plane shape as it impinged on the 
earth's bow shock. This sort of interaction, in the case of the January 11, 
1968 shock, has been studied by Leaning and Chao (1971) and earlier, in 
its general considerations, by Dryer et al. (1967)- The qualitative 
view developed here is in essential agreement with these previous 
investigations . 
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